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Abstract 

We investigate the phonon mechanism for the orientational pinning of Wigner crys- 
tals and stripes in two-dimensional electron layers in GaAs matrices. We find the 
orientation of bilayer Wigner crystals on the (001), (111), (Oil) and (311) interfaces 
versus the interlayer distance and determine the regions of parameters, where poly- 
domain structures can emerge. For the stripe states in electron layers situated close 
to the (001) surface we show that the interference between the piezoelectric and 
deformation potential interaction may be responsible for the preferable orientation 
of the stripes along the [110] direction. For the bilayer system on the (001) interfaces 
we predict the suppression of the resistance anisotropy. 
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1 Introduction 



Due to the number of technological advantages AlGaAs heterostructures are 
the most widely used systems to construct two-dimensional electron layers. 
Low level of disorder in this systems make it possible to realize the states with 
nonuniform spatial distribution of electrons. The formation of these states is 
regulated by the Coulomb interaction, which, under certain conditions, forces 
the transition into the Wigner crystal or the stripe state. Since in cubic crystals 
the Coulomb interaction is the isotropic one, the energy of such states does 
not depend, in the first approximation, on their orientation relative to the 
crystallography axes of the host matrix. Nevertheless, experimental studies of 
the resistance anisotropy in high Landau levels [1-3] (which is considered as 
an indication of the stripe state) show that certain orientational symmetry- 
breaking mechanisms exist. 
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In this paper we consider the electron-phonon interaction as a possible origin 
of the native anisotropy. The main source of the anisotropy is the piezoelec- 
tric interaction, which remains anisotropic in cubic systems. This mechanism 
was studied by Rashba and Sherman [4,5]. In [4,5] the influence of the piezo- 
electricity on the structure and orientation of monolayer Wigner crystals in 
isotropic piezoelectrics was investigated. Since in GaAs crystals the anisotropy 
of the elastic moduli is quite large, the results of [4,5] cannot be applied di- 
rectly to them. In this paper, using the approach lightly different from [4,5], 
we consider the model, where the anisotropy of the elastic moduh of GaAs 
is taken into account. In section 2 the model is applied to the study of ori- 
entation of bilayer Wigner crystals. Technical details of the calculations were 
given in Ref. [6]. Here we mainly concentrate on the results. For the bilayer 
systems on the (001), (111), (Oil) and (311) interfaces we find the stable and 
metastable orientations of the Wigner crystal versus the interlayer separation 
and determine the regions of parameters, where polydomain Wigner crystal 
structures can be realized. 

In section 3 we address the problem of stripe orientation in high Landau levels. 
It was shown in [7] that in an electron layer parallel to the (001) crystallogra- 
phy plane the piezoelectricity causes the preferential orientation of the stripes 
along the [110] or the [110] direction, while the piezoelectric mechanism alone 
does not explain further lowering of the symmetry of the interaction (in mag- 
netic field perpendicular to the layer the [110] orientation is only observed). 
We extend the model of Ref. [7] and take into account the piezoelectric as 
well as the deformation potential interaction. It is shown that for the elec- 
tron layers situated near the surface of the sample these two channel of the 
electron-phonon interaction are strongly interfere and due to this effect the 
C^y symmetry of the electron-electron interaction is reduced to the C2v one. 
We also analyze the orientation of the stripes in the electron layers parallel to 
the (111), (Oil) and (311) crystallography planes, and in bilayer systems. 



2 Piezoelectric mechanism of orientation of bilayer Wigner crystals 



The electrostatic potential ip of an electron in a piezoelectric matrix can be 
found from the solution of the following equations: 

VD = 47re(5(r), (1) 
= 0, 2 



where = —edi^p — '^T^l3i^ki'Uki is the electric displacement vector, — 
KkimU'im~Pi.ikdi'p, the stress tensor, e, the dielectric constant, Xikim, the elastic 
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moduli tensor, (3i^ki, the piezoelectric moduli tensor, Uik, the strain tensor. 

Solving Eqs. (1,2) one can present the electron-electron interaction potential 
in the form: 

V{V) ^^-x-Y: E AnmYnmi^r, <f>r) + O(x') (3) 

^'^ n>0 |m|<n 

where Ynm{(^, 4>) ai'e the spherical harmonics, x = ^m/^Cu ( here and below the 
cubic symmetry of the host matrix is assumed). For the case of an anisotropic 
crystal (cn 7^ Cu + 2C44) the coefficients Anm are computed numerically. 

The second term in r.h.s. of Eq. (3) describes the second order correction to 
the electron-electron interaction caused by the virtual exchange of acoustic 
phonons. This term contains the dependence on the direction of the bond. 

The lattice sum for the interaction (3) can be reduced to the rapidly convergent 
form (see Ref. [6]). Using the rapidly convergent form we have computed 
the energy of the bilayer Wigner crystal and found its stable and metastable 
orientations versus the parameter 77 = d^Jn/2 {d is the interlayer distance, 
and n is the electron density). In principle, the piezoelectric interaction may 
influence on the symmetry of the electron lattice as well, but for GaAs, where 
the parameter x is small 2 • 10~^), this effect is negligible, and the phase 
diagram is almost the same as in case of pure Coulomb interaction [8] . 

For the calculations we use the elastic moduli of GaAs: Cn = 12.3-10^° N/rn^, 
= 5.7 • 10^° N/m?, C44 = 6.0 • 10^° N/rn^. We would like to note that the 
results do not change qualitatively under the 10% variation of Cj^. Therefore, 
one can neglect a small difference of Cik in pure GaAs and in AlGaAs. The 
results of the calculations for the (001), (111), (Oil), and (311) interfaces are 
presented in Fig. 1, where the most preferable orientation of the shortest prim- 
itive lattice vector is shown. One can see that in general case the orientation of 
the bilayer Wigner crystal is sensitive to the interlayer separation. A common 
feature of the dependences presented is a jump-like reorientation under the 

1- st order type transition from the rhombic to the hexagonal phase. For the 
(111) and (Oil) interfaces a step-hke reorientation is also observed under the 

2- nd order type transition from the rectangle to the cubic phase. 

It is instructive to determine the regions, where a polydomain structure of the 
bilayer Wigner crystal may emerge. As it follows from Fig. 1 the polydomain 
structure is expected for the rectangle, rhombic and hexagonal phases in the 
(001) oriented bilayers, for the rectangle, square and rhombic phases in the 
(111) bilayers, for the rectangle phase in the (Oil) bilayers, and for all phases 
(except the one-component hexagonal) in the (311) bilayer. In the last case 
the rectangle phase may show such a behavior at 77 > 0.17, and the cubic 
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Fig. 1. Direction of the shortest in-layer electron-electron bond versus the inter- 
layer separation for the bilayer Wigner crystal on the (001), (111), (Oil), and (311) 
interfaces. Solid curves correspond to the stable orientation; dotted curves, the 
metastable one. The Roman numbers indicate the regions, where different types 
of the electron lattice are realized: I, staggered rectangle; II, staggered cubic; III - 
staggered rhombic; IV, staggered hexagonal (the one-component hexagonal phase 
region is too narrow to be shown in this scale) 



phase - at ?7 > 0.47. 



3 Piezoelectric mechanism of stripe orientation 



In this section we consider the role of the electron-phonon interaction in the 
orientation of stripe states in high Landau levels. We consider an electron 
layer situated close to the (001) surface of the sample and take into account 
two channel (piezoelectric and deformation potential) of the electron-phonon 
interaction. For the (001) surface the matrix elements of the interaction of the 
electrons with surface acoustic modes are in-phase for these two channels (see 
[9]), and the phonon mediated interaction contains the interference term. The 
interference term may lower the symmetry of the electron-electron interaction. 
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The potential energy is given by the expression 



U 



-I 



(fri 
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z>0 



E 



(4) 



where E is the electric field, and the deformation potential interaction is chosen 
in the form 



Udef = j d^rAp{uxx + Uyy)5{z + a) . 



(5) 



Here A is the deformation potential constant, p, the electron density, a, the 
distance between the electron layer and the surface. Since we consider the 
model of a zero-thickness electron layer, the interaction with u^z deformations 
is not included in (5). 

The quantities D and a-i^ satisfy the equations VD = 0, dkU-ik = 0. The bound- 
ary conditions at the free surface are the continuity of Dz and E^, Ey, and 
the vanishing of aiz- At 2; = —a the quantities and are discontinuous: 



-a+O 



z=-a+0 



-a-0 



z=—a—0 



(6) 
(7) 



where F = A{dxP, dyp, 0) is the tangential force applied to the unit area of the 
interface. This force is induced by the deformation potential interaction. We 
should note that discontinuity of the stresses is the consequence of the zero- 
thickness approximation for the electron layer. Using the boundary conditions 
specified one can reduce the energy (4) to the form 

U = lj dh{ep^ - UiFi) , (8) 

where the electrostatic potential (/? and the displacement field u are taken at 
z = —a. The energy (8) can be presented as a sum of a pure Coulomb part 
(which is isotropic) and a phonon part C/ph (which includes the isotropic as 
well as the anisotropic terms). 

For simplicity, we consider the unimodal approximation for the electron den- 
sity modulation p — po cos qr^;, where the absolute value of q is determined by 
the period I of the stripe structure {q — 27r/Z) and its direction is perpendicular 
to the axis along which the stripes are aligned. 

We concentrate on the anisotropic GaAs crystal and compute the values of 
(f and Ui at z — —a numerically. We use the parameters eu — 0.15 C/m^, 



5 



[100] [110] [010] [110] [100] 



Fig. 2. Ptionon contribution to the energy of ttie stripe phase (in units of xUc) near 
the (001) surface versus the direction of the stripes at a/l =1.0, 0.6, 0.5, 0.4, 0.2 
(from top to bottom). 

A = 7.4 eV, e = 12.5, 1 = 2- 10^ A, and the elastic moduli given above 
for the calculation. For the (001) surface the dependence of the energy f/ph 
on the direction of stripes is shown in Fig. 2. The values of f/ph are given in 
units of xUc (per unit area), where Uc = 'Ke'^pl/2eq is the Coulomb energy 
at a/Z » 1 Our calculations show that at a/Z > 0.23 the global minima are 
reached for two directions symmetrically deviated from the [110] axis on a 
small angle. In this case one can expect that a domain structure is formed, 
and, in average, the system should demonstrate the minimum resistance in the 
[110] direction and the maximum resistance in the perpendicular direction. 
The calculations predict the largest resistance anisotropy at a/l — 0.4. If 
there is no other sources for the native anisotropy, the high and low resistance 
directions may alternate at a/l ^ 0.23. At a/l > 1 the interference term 
becomes exponentially small and the 4-fold symmetry is restored. 

For po of order of the average density of electrons at the high half-filled Landau 
level and a/l = 0.4 we estimate the absolute value of the anisotropy (the 
energy difference of the [110] and [110] oriented stripes) as 0.7 mK per electron, 
which is comparable with the values given by the other possible mechanism 
of the anisotropy [10]. The mechanism [10] is connected with the effective 
mass anisotropy caused by the asymmetry of potential confining the electron 
in the quantum well. From our estimate we conclude that both mechanisms 
may work in parallel and in case of the symmetric confining potential [3] the 
native anisotropy can be accounted for the phonon mechanism. Note that the 
effect considered depends significantly on the distance between the electron 
layer and the surface. 

Assuming the phonon mechanism plays am important role in the orientational 
pinning of the stripes, it is interesting to analyze the orientation in electron 
layers parallel to some other crystallography planes. The results of the cal- 
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Fig. 3. Stripe energy anisotropy for tlie (111), (Oil) and (311) surfaces. Solid curve 
- a/l = 1, daslied curve - a/l = 0.4, dotted curve - a/l = 0.2. ?7ph is in x^c units. 

culations for the (111), (Oil) and (311) plane are shown in Fig. 3. In these 
cases, as it follows from the results presented, the orientation of the stripes 
does not depend on the ratio a/l. For the (311) surface the phonon mecha- 
nism predicts the [233] orientation. This prediction is in agreement with the 
experiment [11]. The absolute value of the anisotropy for the (311) surface is 
in several times large then for the (001) one. For the (111) plane the global 
minima correspond to three different directions of the stripes ([21 1], [1 2 1] and 
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Fig. 4. Stripe energy anisotropy (in x^c per one layer) for the bilayer system near 
ttie (001) surface at a/l=0.4. Solid curve - d/l = 0.5, dashed curve - d/l = 1, dotted 
curve - d/l = 2 

[112]). In tills case the polydomain structures may not sliow any resistance 

anisotropy. For the (Oil) plane the domains with [111] and [111] orientation 
are energetically preferable, and the polydomain structure will show the re- 
sistance anisotropy with the [Oil] low resistance direction. The anisotropy is 
expected to be small due to a large angle (~ 70°) between two preferable 
orientations. 

To complete our study we consider the orientation of stripes in the bilayer 
systems. Let there are two electron layers on the distances a and a-\-d from the 
surface, and the electron density modulation in the one layer is pi = po cos qVpi, 
and in the other layer is p2 — — pocosqrpj. For the (001) surface the change 
of the anisotropy under variation of the ratio d/l is shown in Fig. 4. On can 
see that in the bilayer systems with a polydomain structure the angle between 
the directions of stripes in different domains becomes larger for smaller d/l. 
At d/l = 0.5 these directions are almost perpendicular to each other. It means 
that the resistance anisotropy becomes very small. For the (111), (Oil) and 
(311) crystallography planes we do not find any significant changes in the 
stripe orientation under the variation of the parameter d/l. In these cases the 
bilayer and monolayer systems will demonstrate the same orientation. 

To conclude this section, we would like to point out on the questions, which, 
in our opinion, require further experimental investigation. 1. Is there a depen- 
dence of the resistance anisotropy on the distance between the electron layer 
and the surface? 2. Do the electron layers on other interfaces demonstrate 
the resistance anisotropy? 3. Is this effect modified in bilayer systems? Ex- 
perimental study of these question may be important for establishing of the 
mechanism for the orientational symmetry breaking in high Landau levels. 
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